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ABSTRACT 

The effect of an arbitrary change of frame on the structure of turbulence 
models is examined from a fundamental theoretical standpoint. It is proven, 
as a rigorous consequence of the Navier-Stokes equations, that turbulence 
models must be form invariant under arbitrary translational accelerations of 
the reference frame and should only be affected by rotations through the in- 
trinsic mean vorticity. A direct application of this invariance property 
along with the Taylor-Proudman Theorem, material frame-indifference in the 
limit of two-dimensional turbulence and Rapid Distortion Theory is shown to 
yield powerful constraints on the allowable form of turbulence models. Most 
of the commonly used turbulence models are demonstrated to be in serious vio- 
lation of these constraints and consequently are inconsistent with the Navier- 
Stokes equations in non-inertial frames. Alternative models with improved 
non-inertial properties are developed and some simple applications to rotating 
turbulent flows are considered. 


This research was supported by the National Aeronautics and Space Admini- 
stration under NASA Contract No. NAS1-18107 while the author was in residence 
at the Institute for Computer Applications in Science and Engineering (ICASE), 
NASA Langley Research Center, Hampton, VA 23665. Partial support was also 
provided by the Office of Naval Research under Contract No. N00014-85-K-0238. 


i 


1. INTRODUCTION 


Turbulence plays a fundamental role in a variety of physical systems 
which evolve in non-inertial frames of reference. Various types of fluid 
machinery and geophysical systems (e.g., gas turbines, propellers, ocean cur- 
rents, and atmospheric weather fronts which can have a profound effect on our 
daily lives) generate turbulence in non-inertial reference frames that are 
undergoing time-dependent rotations and translations relative to an inertial 
framing. Consequently, it is essential that a clear understanding of such 
non-inertial effects on turbulence be gained if these physical systems are to 
be modeled properly. Furthermore, due to the analogy between rotations and 
curvature, a physical model which does not properly account for non-inertial 
effects is likely to yield erroneous predictions for problems involving curva- 
ture in inertial frames of reference. 

To date, there have been no comprehensive studies of non-inertial effects 
on turbulence modeling based on a rigorous analysis of the Navier-Stokes equa- 
tions. Most of the previous studies consisted of rigorous mathematical 
analyses of the highly simplified limiting case of two-dimensional turbulence 
(see Speziale 1981, 1983) or more applied studies of three-dimensional turbu- 
lence where the effects of rotations of the reference frame were accounted for 
by a variety of ad hoc empiricisms (c.f., Majumdar, Pratap, and Spalding 1977, 
Howard, Patankar, Bordynuik 1980, and Galmes and Lakshminarayana 1983). There 
have been several studies applying second-order closure models to turbulent 
flows in rotating frames which are substantially less empirical in nature 
(c.f., Mellor and Yamada 1974, So 1975, So and Peskin 1980, and Launder, 
Tselepidakis , and Younis 1987). However, it was recently proven by Speziale 
(1985) that these particular second-order closure models are fundamentally 
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inconsistent with the Navier-Stokes equations in a rapidly rotating frame. 
Consequently, such models cannot be applied to turbulent flows in arbitrary 
non-inertial frames of reference without the need for making ad hoc adjust- 
ments in the constants. Since direct numerical simulations of such turbulent 
flows, at the high Reynolds numbers and in the complex geometries of scien- 
tific and engineering interest, will not be possible for at least the next 
several decades, it is essential that turbulence models be developed whose 
properties in non-inertial frames of reference are consistent with the Navier- 
Stokes equations. This forms the raison d'etre of the present study. 

In this paper, it will be proven that turbulence models should be form 
invariant under arbitrary translational accelerations of the reference frame 
relative to an inertial framing (i.e., the exact invariance group of turbu- 
lence models is the extended Galilean group). Rotations of the reference 
frame will be shown to affect turbulence models only through the intrinsic 
mean vorticity. It will be shown that these rotationally-dependent non- 
inertial effects must vanish for a two-dimensional turbulence (i.e., material 
frame-indifference in the limit of two-dimensional turbulence; see Speziale 
1981, 1983) and should be consistent with Rapid Distortion Theory and the 
Taylor-Proudman Theorem (c.f., Greenspan 1968). A systematic application of 
these ideas will be shown to provide powerful constraints on the allowable 
form of turbulence models. A variety of the popularly used turbulence models 
(e.g., zero, one, or two equation turbulence models along with second-order 
closures) will be shown to be in serious violation of these constraints which 
can give rise to spurious physical results in rotating frames. Improved two- 
equation models and second-order closure models will be presented along with 
some brief applications to rotating turbulent flows. 
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2. CONSTRAINTS ON TURBULENCE MODELS IN NON— INERTIAL REFERENCE FRAMES 

We will consider the incompressible turbulent flow of a homogeneous 
viscous fluid governed by the Navier-Stokes and continuity equations which 
take the form (c.f., Batchelor 1967) 

+ v • Vv = -VP + vV^v - ftxx-ax(axx)-ft^-20xv (1) 

V • v = 0 (2) 


in an arbitrary non-inertial reference frame (see Figure 1). In Equations 
( 1 )— (2 ) , v is the velocity vector, P is the modified pressure, ft(t) is the 

rotation rate of the non-inertial frame relative to an inertial framing, 

• 

UgCt) is the translational acceleration of the origin of the non-inertial 
frame relative to an inertial framing, and v is the kinematic viscosity of 
the fluid. It should be noted that the Navier-Stokes equations are altered by 
the presence of four frame-dependent terras on the right-hand-side of (1) 

which, respectively, are referred to as the Eulerian, centrifugal, transla- 
tional, and Coriolis accelerations. The continuity equation is frame- 

indifferent, i.e., it has no non-inertial terms and, hence, is of the same form 
in all frames of reference independent of whether or not they are inertial. 

As in the usual treatments of turbulence, the velocity field v and 

pressure P will be decomposed into ensemble mean and fluctuating parts as 

follows : 

v = v + u, P = P‘+p (3) 


where 


i N 

v ■ iim l v 

N-k» a=l 


(a) 


1 


N 


7 = lim ^ l P 
N-h» a=l 


(a) 


» 


(4) 
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are ensemble averages taken, in practice, over a large number of N realiza- 
tions of the turbulence (c.f., Hinze 1975). For a statistically steady or 
homogeneous turbulence, the ergodic hypothesis can be invoked and time 
averages or spatial averages, respectively, can be substituted. The mean 
velocity v and mean pressure P* are solutions of the Reynolds equation 
and mean continuity equation which, respectively, take the form 

+ v • Vv = -VP + vV 2 v + V*T-axx-ftx(Qxx)-ij - 211 x v (5) 
9 1 o 

V • v = 0 (6) 

in any arbitrary non-inertial reference frame where 

T “ -5u (7) 

is the Reynolds stress tensor. Equations ( 5 )— (6 ) are obtained by substituting 
the decomposition (3) into the Navier-Stokes equations and then taking an en- 
semble average. The fluctuating velocity u and fluctuating pressure p are 
solutions of the following equations (valid in an arbitrary non-inertial 
frame) : 

1^- +▼ • Vu = -o • Vu - u • Vv - Vp + vV 2 u - ? • t - 28 x a (8) 

9 1 

V • u * 0 (9) 

which are referred to as the fluctuating momentum and fluctuating continuity 
equation, respectively. Equations (8 )— (9 ) are obtained by subtracting 
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Equations (5 )— (6 ) from Equations (l)-(2), respectively. From Equation (8), it 
is clear that the evolution of the fluctuating velocity u (for a given mean 
velocity field) is only affected by the reference frame through the Coriolis 
acceleration 20 x u. Eulerian, centrifugal, and translational accelera- 
tions only have an indirect effect on the fluctuating velocity through the 
changes that they induce in the mean velocity. 


At this point, the concepts of the Oldroyd derivative and intrinsic 
vorticity will be introduced. The Oldroyd derivative of the fluctuating 


velocity is defined by 


or 5 §? + " 


and represents the frame-indifferent convected time rate of u following the 
mean velocity with respect to both position and orientation. Unlike the sub- 
stantial derivative Du/Dt = 3u/3t + v • Vu, the Oldyroyd derivative is 
independent of the observer; relative to any two independent non-inertial 
reference frames x and x (whose motions can differ by an arbitrary time- 
dependent rotation and translation) the Oldroyd derivative of a given fluctu- 
ating velocity field is the same, i.e., 

Du* Du 

-V - — . do 

Dt Dt 

The intrinsic vorticity W (also referred to as the absolute or potential 
vorticity) is the vorticity relative to an inertial framing and is obtained by 
adding 2Q to the local vorticity ta = V x v in the non-inertial 


frame 


Hence, we have (c.f., Tritton 1977) 
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W = la + 2Q . 


( 12 ) 


By taking the ensemble mean and dual of (12), the intrinsic mean vorticity 
vector and tensor are obtained which, respectively, are given in component 
form by the equations 

\ = \ + 2 \ (13) 



k i 


m£k m 


(14) 


where uT^ = (3"v^/3x^ - 3v^/3x^) is the local mean vorticity tensor and 

e m£k P ermutat ^ on tensor. A straightforward substitution of (10) and 

(14) into (8) gives rise to the alternative form of the fluctuating momentum 
equation 

Du _ 2 

— — = -u»Vu-2(s + W)*u-Vp + vVu-V»T (15) 


which is valid in any non-inertial frame of reference. In (15), W is the 
intrinsic mean vorticity tensor and "5 is the mean rate of strain tensor 
whose components are given by 


kfc 


1 3v k 9v £ 

7 ( JZ7 + 

i k 


(16) 


in all frames of reference independent of whether or not they are inertial. 
As a result of (15), it is clear that the evolution of the fluctuating 
velocity (and higher-order moments constructed from it) only depend on the 
reference frame through the intrinsic mean vorticity W. 
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The Reynolds and continuity equations (5)-(6) are not closed as a result 
of the additional unknowns represented by the six components of the Reynolds 
stress tensor t^. Tn virtually all previous studies of turbulence model- 
ing beginning with Boussinesq, it was tacitly assumed that the Reynolds stress 
tensor is uniquely determined by the global history of the mean velocity 
field. This assumption is generally consistent with the Navier-Stokes equa- 
tions in an inertial framing as pointed out by Lumley (1970). Hence, in an 
inertial frame of reference, we have 

T = x[v(x',t'); x,t] x^D, t' € (-<», t) (17) 


where D is the fluid domain and a bracket [•] denotes a functional (i.e., 
any quantity determined by the global history of a function). It should be 
noted at this point that both t and v are kinematical quantities whose 
transformation properties under a change of frame are mathematically deter- 
mined. To be more specific, given that x is an inertial frame and x* is 
an arbitrary non-inertial frame, it is a simple matter to show from basic 
kinematics that (c.f., Speziale 1979) 

* 

T = T (18) 


— * _ 

v = v - Q x x - Uq. 


(19) 


Consequently, once the inertial form of (17) is specified its non-inertial 
form is automatically determined . It thus follows that if the non-inertial 
form of (17) is incorrect, its inertial form must also be incorrect since the 


two are not independent 
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Since x is uniquely determined from the fluctuating velocity (i.e., 
x = - uu) , it is clear that its invariance group must be at least as large as 
the invariance group of the fluctuating momentum equation (8) (of course, (17) 
constitutes a special solution of (8) and, hence, it could have a larger in- 
variance group; c.f., Rosen 1980). Consequently, turbulent closure models 
for x must be at least form invariant under arbitrary translational ac- 
celerations of the reference frame. Hence, Equation (17) must transform in 
the form invariant manner 

xfv*(x',0; x ,t] = x[v(x',t'); x,t] , x'€D, t'e(-«,t) (20) 

under the extended Galilean group of transformations 

x* * x + c(t) (21) 

• •• 

where Uq = -c is the translational acceleration of the non-inertial frame 
relative to an inertial framing. Constraint (20) would, for example, forbid 
turbulent closure models from having any explicit dependence on the mean 
acceleration Div/Dt . Since any dependence on the rotation rate of the 
reference frame must arise from the intrinsic mean vorticity, it follows that 
in an arbitrary non-inertial frame , turbulent closure models for x must be 
of the general form 

T - x[v(x',t'), W(x',t'); x,t] x'€D, t' €(-«,t). (22) 

Here it is understood that the explicit functional dependence on v in (22) 
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is frame-indifferent (i.e., does not contain any terms which depend on the 
motion of the reference frame relative to an inertial framing). Frame- 
dependence only enters implicitly through the intrinsic mean vorticity W. 
Examples of one-point fields that are frame-indifferent functionals of the 
mean velocity include "S and D^/Dt; two-point fields that are frame- 

indifferent include the vorticity difference 5(x',t) - S(x,t) and its 
Oldroyd derivative. 

Although three-dimensional turbulent closure models can be frame- 
dependent through the intrinsic mean vorticity tensor W, it has been shown 
recently that such models must become frame-indifferent in the limit of two- 
dimensional turbulence (see Speziale 1981, 1983). By a two-dimensional turbu- 
lence we mean a turbulent flow where the fluctuating velocity u is of the 
plane two-dimensional form 

u = u (x,y,t) i + u (x,y,t) j . (23) 

x y 

Consistent with this two-dimensional assumption, the angular velocity of the 
reference frame must be of the form Q = fik so that the mean velocity in- 
duced by it is comparably two-dimensional. For such a flow, the Coriolis 
, acceleration in the fluctuating momentum Equation (8) is derivable from a 

scalar potential as follows (see Speziale 1981, 1983) 

2fl x u = V(2fli|>) (24) 

since, as a general solution of the two-dimensional continuity equation, the 
fluctuating velocity can be written in the stream function form 
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1 u = Vip x k. (25) 

Consequently, the Coriolis acceleration (which constitutes the only non- 
inertial effect in Equation (8)) can be absorbed into the fluctuating pressure 
in a two-dimensional turbulence leaving the fluctuating velocity unaffected. 
Consistent with this result, the dependence of the Reynolds stress tensor on 
the intrinsic mean vorticity (which characterizes these Coriolis effects) must 
vanish in the limit of two-dimensional turbulence rendering the model frame- 
indifferent. Thus, in the two-dimensional limit, turbulence closure models 
for the Reynolds stress tensor must be of the same form 

x « T[v(x',t'); x,t] x'£D, t'€(-®,t) (26) 

independent of whether or not the reference frame is inertial. This invari- 
ance property is referred to as the principle of material frame-indifference 

i 

I 

In the limit of two-dimensional turbulence (see Speziale 1981, 1983). 

The limit of two-dimensional turbulence constitutes a real physical limit 
which can be approached by any statistically steady turbulence, sufficiently 
far from solid boundaries, in a rapidly rotating framework (a direct conse- 
quence of the Taylor-Proudman Theorem; c.f., Tritton 1977). The Taylor- 
Proudman Theorem in its classical form states that steady inviscid flows in a 
rapidly rotating framework are two-dimensional, i.e., are independent of the 
coordinate along the axis of rotation of the fluid. Of course, the Taylor- 
Proudman Theorem holds in an excellent approximate sense for most laminar vis- 
cous flows provided that the flow is sufficiently far removed from solid 
boundaries where Ekman layers can develop. For a statistically steady 

i 

I 
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turbulent flow In a rotating frame, the Reynolds equation takes the form 

v • Vv = -VP + vV^v M • i - 2D x v (27) 

where the centrifugal acceleration has been absorbed into the modified pres- 

sure F. The associated mean vorticity transport equation, obtained by 
taking the curl of Equation (27), is given by 

v • Vu» = u • Vv + vV^u + V x (V • t) + 2fl • Vv (28) 

where w = V x v is the local mean vorticity in the rotating frame. If we 

let Q = Qk, Equation (28) can be written in the alternative form 

_ 1 _ _ _ _ 0 — 

k • Vv ■ (v • Vu - # • Vv - v7 u - 7 x V « t). (29) 

In the limit as &•*•», Equation (29) reduces to 

k • V v e = 0 (30) 

a Z 

for a statistically steady turbulence. Sufficiently far from solid 

boundaries, Equation (30) has the simple solution 

v = v(x,y) (31) 


and, thus, the mean velocity field for a statistically steady turbulence in a 
rapidly rotating frame must be two-dimensional. If the flow is confined by 
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boundaries normal to the axis of rotation (at distances sufficiently far re- 
moved from the flow region being considered), the mean velocity field will be 

of the two component form v = v (x,y)i + v (x,y)j (see Tritton 1977). The 

x y 

same type of two-dimensionalization will hold for the Reynolds stress tensor 
in an approximate sense since the filtered velocity satisfies an equation of 
the form of (28) and the large scales of turbulence contain most of the 
energy. This is the turbulent generalization of the classical Taylor-Proudman 
Theorem which has been verified experimentally to hold in an excellent approx- 
imate sense provided that the Rossby number Ro = Vq/Q£q << 1 (where 
and vq are the length and velocity scales of the turbulent fluctuations)* 
Any turbulence model which does not yield such a Taylor-Proudman reorganiza- 
tion in a rapidly rotating framework is fundamentally inconsistent with the 
non-inertial form of the Navier-Stokes equations. 

The last constraint that will be considered in this section involves the 
application of Rapid Distortion Theory (RDT) to turbulence suddenly subjected 
to a strong rotation. Since the Taylor-Proudman Theorem serves primarily as a 
constraint on the large energy containing eddies in a rapidly rotating frame, 
RDT will be used as a constraint on the small scale turbulence which is not 
too far removed from isotropy. Hence, RDT for an initially isotropic turbu- 
lence which is suddenly subjected to a rapid rotation will be considered. For 
this problem, the initial Reynolds stress tensor and dissipation rate tensor 
are of the isotropic form 


T 


ij 


7 Vij 


(32) 



2 

1 


e 0 6 iJ 


(33) 
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respectively, where Kq is the initial turbulent kinetic energy and is 
the initial dissipation rate of the turbulence (it should be noted that 

= 2 v (3u^/3x m )(3Uj /<* x m ) an< i e = y Rapid Distortion Theory 
predicts that the turbulence undergoes an Isotropic linearly viscous decay 
(c.f., Reynolds 1987). More specifically, according to RDT, the Reynolds 
stress tensor and dissipation rate tensor are of the isotropic form 

x ij --f K(t) V D ij - ! e<t)S u (34) 

at a later time t > 0 after the turbulence has been subjected to the rapid 
rotation. Here, the turbulent kinetic energy and dissipation rate are deter- 
mined from the energy spectrum E(k,t) as follows (c.f., Hinze 1975 ) 


00 


K(t) = / 
0 

E(k,t)dk 

(35) 

00 

e(t) = 2v / 

0 

k 2 E(k,t)dk 

(36) 


which are valid for an isotropic turbulence. Rapid Distortion Theory predicts 
that the energy spectrum undergoes a linearly viscous decay, and thus at any 
later time t > 0: 

E(k,t) = E(k,0)exp(-2vk 2 t) (37) 

(it should be noted that In the limit of infinite Reynolds numbers, the energy 
spectrum remains unchanged for finite times t > 0; c.f., Reynolds 1987). Of 
course, RDT is only formally valid for short elapsed times t « K^/e^. 
However, since a rapid rotation destroys the phase coherence needed to cascade 
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energy from the large scales to the small scales (represented by the energy 
transfer term), it would appear that RDT could hold for much larger elapsed 
times for the case of a rapidly rotating isotropic turbulence. This was veri- 
fied by the results of direct numerical simulations of the Navier-Stokes equa- 
tions for isotropic turbulence subjected to a rapid rotation (see Speziale, 
Mansour, and Rogallo 1987). These numerical simulations indicated that the 
rapid rotation suppresses the energy transfer for long time intervals yielding 
results in excellent approximate agreement with the RDT solution specified by 
Equations (34) and (37) for the primary period of the decay (i.e., up to and 
beyond the point where the turbulent kinetic energy has decayed to 10% of its 
initial value). The time evolution of the energy spectrum obtained from a 
128^ direct numerical simulation is shown in Figures 2(a)-(b) for a Reynolds 
number Re^ = 51 and a Rossby number Ro^ = 0.07 based on the initial 
turbulent kinetic energy and Taylor microscale. These computed energy spectra 
were illustrative of a linearly viscous decay during the entire period of the 
computation (i.e., for 0.1 < K/K q < 1.0). The L 2 norm of the anisotropy 
tensor was extremely small (llbll^ < 0.01) for the entire duration of the 
computation and, hence, the rotation had no discernible effect on the isotropy 
of x. These results demonstrate strong agreement with RDT for large 
elapsed times. Numerical results for the decay of the turbulent kinetic 
energy (shown in Figure 3) illustrate that the rapid rotation gives rise to a 
dramatically reduced turbulence dissipation rate due to the disruption of the 
energy transfer from large scales to small scales. It is the opinion of the 
author that these fundamental results (which are important since they capture 
the essential physical features of the reaction of small scale turbulence to a 
rapid rotation for long as well as short elapsed times) should serve as a 
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baslc constraint on turbulence models in rotating frames. Models that are in 
serious violation of these RDT results are likely to give rise to spurious 
physical results in rotating turbulent flows. 
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3. INCONSISTENCY OF EXISTING TURBULENCE MODELS 

As demonstrated in the previous section, the Navier-Stokes equations in a 
non-inertial reference frame place the following basic constraints on the 
allowable form of turbulence models: 

(i) Reynolds stress models must be form invariant under arbitrary trans- 
lational accelerations of the reference frame and should only be affected 
by rotations of the reference frame through the intrinsic mean vorticity. 

(ii) All frame-dependent effects (and thus any dependence on the in- 
trinsic mean vorticity) must vanish in the limit of two-dimensional 
turbulence — a constraint appropriately named material frame-indifference 
in the limit of two-dimensional turbulence. 

(iii) Reynolds stress models must be consistent with the Taylor-Proudman 
Theorem for turbulent flows. This requires that a statistically steady 
turbulence in a rapidly rotating frame (sufficiently far from solid 
boundaries) be two-dimensional . 

(iv) Turbulence models should be consistent with the results of RDT for 
an initially isotropic turbulence subjected to a rapid rotation. This re- 
quires Reynolds stress models to predict that an initially isotropic 
turbulence undergoes an isotropic linearly viscous decay in a rapidly 
rotating frame yielding a substantially reduced dissipation rate. 



-17- 


First, we will note the inconsistency of the more empirical Coriolis 
modified turbulence models with these constraints. These more empirical 
models are characterized by the introduction of coefficients that depend ex- 
plicitly on the rotation rate of the reference frame. For example, in the 
model of Howard, Patankar, and Bordynuik (1980), empirical coefficients in the 
transport equations for the turbulent kinetic energy and dissipation rate in 
rotating duct flow were allowed to depend on the Richardson number 


Ri = 


~2ft (fj - 20) 

(9lT) 2 


(38) 


where ft is the rotation rate of the duct,' u is the mean velocity along 
the axis of the duct, and y is the transverse coordinate. A comparable em- 
pirical model, based on the nonlinear algebraic model of Rodi (1976), was pro- 
posed recently by Galmes and Lakshminarayana (1983) where an implicit depen- 
dence on the Richardson number (given by Equation (38)) was introduced into 
the Reynolds stresses. Such empirical models (which have also been proposed 
by other authors) violate constraint (i) and are thus inconsistent with the 
Navier-Stokes equations. More specifically, rather than an explicit depen- 
dence on the rotation rate ft there should be an implicit dependence on ft 
only through the intrinsic mean vorticity (i.e., the quantity -3ii/9y + 2ft 
for the rotating channel flow under discussion). The recent large-eddy simu- 


lations of Bardina, Ferziger, and Reynolds (1983) for rotating homogeneous 
shear flow demonstrated that the turbulent Reynolds stresses do not scale with 
the Richardson number. 

Eddy viscosity models form the foundation for most of the turbulence 
models that are used by scientists and engineers. These models are of the 
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general form 


where 


D T k£ 2v T S kJt 


D T U “ T kH “ 3 T nmi 6 kl 


(39) 

(40) 


is the deviatoric part of the Reynolds stress tensor and is the eddy 
viscosity in its kinematic form. Equation (39) encompasses a wide variety of 
turbulence models which, by far, are the most commonly used models for the 
solution of practical problems. We will now examine the consistency of a 
variety of popular eddy viscosity models with constraints (i)-(iv) for ro- 
tating turbulent flows. The simplest eddy viscosity models are the zero equa- 
tion models where the turbulent time scale is constructed from the mean velo- 
city gradients and the turbulent length scale is specified algebraically. Two 
such popular models are the Smagorinsky (1963) Model given by 


v = £ 2 (2S S ) 1/2 
T mn in 


(41) 


and the Baldwin-Lomax Model (or vortlcity model) given by 

v = l 2 (w w ) l ^ 2 

i mm 


(42) 


where l is the turbulent length scale which is usually specified empiri- 
cally. The Smagorinsky Model, which is the tensorial generalization of 
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Prandtl's mixing length theory*, arose out of the Geophysical Fluid Dynamics 
community and (as a subgrid scale model) has served as the cornerstone for 
large-eddy simulations. Since v^, only depends on £5, it is frame- 

indifferent for all mean flows and, as such, automatically satisfies con- 
straints (ii)-(iii). However, since it is frame-indifferent in three- 
dimensions as well as in two-dimensions, the Smagorinsky Model is fundamental- 
ly incapable of describing the effects of rotation in retarding the energy 
transfer process (as described in constraint (iv)) which ultimately has an 
effect on pT. However, such effects are primarily manifested in the large 
scales and, consequently, the Smagorinsky model would be satisfactory as a 
subgrid-scale stress model despite the fact that it has undesirable rotational 
properties as a Reynolds stress model (see Bardina, Ferziger, and Reynolds 

1983 and Speziale 1985). 

On the other hand, it will now be shown that the Baldwin-Lomax Model is 
more seriously inconsistent with the Navier-Stokes equations in a rotating 
frame. It should first be noted that the eddy viscosity (42) is specified for 
an inertial framing. However, as alluded to earlier, it follows from basic 
kinematics that (see Speziale 1979) 

* _ _* _ 

x = t, S = S, o> = u - 2Q (43) 

where the starred quantities are relative to an arbitrary non-inertial 

£ 

reference frame x . Hence, given that (42) is the inertial form of the 
Baldwin-Lomax Model, it follows that its non-inertial form is given by 


For a unidirectional turbulent shear flow (with mean velocity v - u(y)i) 
Equation (41) reduces to * it^|du/dy|. 
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v_ = l 2 ( W W ) 1/2 . (44) 

T mm 

While this model is consistent with constraint (i) (i.e., frame-dependent 

effects only enter in through the intrinsic mean vorticity), it is in serious 

violation of constraints (ii)-(iv). To be more specific, in the limit of two- 

2 

dimensional turbulence, W W = (io + 2ft) survives and hence there is a viola- 

m m 

tion of material frame-indifference in the limit of two-dimensional turbulence 
(i.e., in the two-dimensional limit, any dependence on 0 must vanish for 
there to be consistency with the Navier-Stokes equations). Furthermore, since 
according to (44), v T -► 00 as ft -*■ «, the Baldwin-Lomax Model predicts that 
there is an increase in turbulent dissipation corresponding to an increase in 
the rotation rate of the framing which violates constraint (iv). This un- 
bounded growth of v^, as ft -*• also gives rise to the violation of 

r\ 

constraint (iii) — the Taylor-Proudman Theorem. For large ft, ■ 2$rft and 
hence in the limit as ft -*• «, (29) reduces to 

k • Vv = -2V x [V • (JL 2 ?)] (45) 

with the implication that 3v/3z is not necessarily zero (i.e., v * v(x,y) 
for any statistically steady turbulent flow sufficiently far from solid 
boundaries) in violation of the Taylor-Proudman Theorem. It is thus clear 
that vorticity models such as the Baldwin-Lomax Model are likely to yield un- 
physical results for turbulent flows involving strong rotational strains and, 
consequently, do not form a general foundation for either a Reynolds stress or 
subgrid scale stress model. 



- 21 - 


One-equation models involve the solution of an additional transport equa- 
tion for the turbulent kinetic energy. The eddy viscosity for such models is 
of the form 


v 


T 


K 1/2 Jt 


(46) 


where K is the turbulent kinetic energy (obtained from its modeled transport 
equation) and i is an appropriate length scale of turbulence which is 
usually specified empirically based upon the particular flow geometry under 
consideration (see Cebeci and Smith 1974 and Rodi 1984 for a survey of such 
models). Since the transport equation for K is frame-indifferent and 
since i is usually only specified based on the geometry of the flow con- 
figuration, such models are identically frame-indifferent. Due to the fact 
that they satisfy material frame-indifference in three-dimensional turbulent 
flows as well as in two-dimensional turbulent flows, they are unable to pre- 
dict the reduction in turbulence dissipation that results from the application 
of a strong rotation (i.e., such one-equation models are generally consistent 
with constraints (i)-(iii) but in serious violation of constraint (iv)). The 
same precise criticism can be leveled against two-equation turbulence models 
among which the K-e model has become extremely popular during the past 
decade. In the K-e model, the eddy viscosity is represented by 



where C * 0.09 is an empirical constant, K is the turbulent kinetic 
energy, and e is the turbulent dissipation rate. In the K-e model, K 
and e are determined from modeled versions of their transport equations 
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which are usually of the form (see Kanjalic and Launder 1972) 
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(49) 


where D/Dt = + v • V and Cj - C 4 are empirical constants . Equations 

(48) - ( 49 ) are of the same form independent of whether or not the reference 
frame is inertial. Consequently, the K-e model is frame-indifferent for all 
flows thus making it impossible for this model to account for the reduction in 
dissipation that occurs in rotating isotropic turbulence as well as in other 
rotating turbulent flows (i.e., the model is in violation of constraint 
(iv)). Furthermore, the inability of the K-e model to accurately predict 
normal Reynolds stress differences in turbulent flows of engineering 
importance (see Speziale 1987) can be exacerbated further in rotating flows 
where Coriolis effects usually give rise to stronger such anisotropies. 

Problems of a similar nature exist with second-order closure models. All 
of the commonly used second-order closure models are of the general form (see 
Speziale 1985) 
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is the frame-indifferent Oldroyd derivative of x and is a constant 

(which arises from the rotationally dependent part of the rapid pressure- 
strain correlation). Here, is a function of the variables x, Vx, 

and J l which arises from the modeling of the triple velocity and pressure- 

diffusion correlations whereas ir^ is a function of x, S, and J l which 
arises from the irrotational part of the production terms and the modeling of 
the slow pressure-strain and dissipation rate correlations. This general form 
(50) encompasses the Launder, Reece, and Rodi (1975) model, the models of 
Lumley (1978), and the Rotta-Kolmogorov model (see Mellor and Herring 1973). 
In the former two models, the length scale of turbulence i is taken to be 
of the form 

Z = K 3/2 /e (52) 


where the dissipation rate e is determined from a modeled transport equa- 
tion which is of the same general frame-indifferent form as (49). Analogous- 
ly, in the Rotta-Kolmogorov model the length scale l is obtained from the 
transport equation (see Mellor and Herring 1973): 


L t(v + Bi k 1/2 i> ji- <«) + 6 2 k 3/2 « H-] 


_ _ i 0 v 3/2 

+ B 3* T ij -5x7 ' 6 4 K 

J 


(where 0 ^ &j> an d 84 are empirical constants) which is of the same 
form in all frames of reference independent of whether or not they are iner- 
tial. Consequently, since D^j^/Dt, and v are f rame_ 
indifferent along with the transport equations for e and Z, it follows 
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that frame-dependence in the commonly used second-order closures arises 
exclusively from the term 


(C. - 2)(x, W. 

1 km Jim 


+ x W ). 

Am km 


(54) 


Since is a constant that does not equal 2 (in the Launder, Reese, and 

Rodi model, Cj = 0.6 whereas in the Rotta-Kolraogorov model, = 0) and 

since x, W. + t„ W. does not generally vanish in a two-dimensional 

km Am Am km 6 J 

turbulence, it follows that material frame-indifference in the limit of two- 
dimensional turbulence is violated. This inconsistent dependence of (50) on 
W also gives rise to a violation of the Taylor-Proudman Theorem in problems 
of engineering and geophysical interest since the constraint 


Aim 


1 


V x (V 


x) = 0 


(55) 


(which is a necessary condition for the Taylor-Proudman Theorem) is violated 
in statistically steady turbulent flows by these second-order closures. 
Furthermore, since the transport equations for e and A are frame- 
indifferent in the commonly used second-order closures, they are unable to 
account for the reduction in dissipation (and the associated change in length 
scales) in rotating isotropic turbulence. Thus, for turbulent flows in a 
rapidly rotating frame, the commonly used second-order closure models are in 
rather serious violation of the Navier-Stokes equations. Although Launder, 
Tselepidakis , and Younis (1987) were able to get reasonable correlation with 
experiments on rotating channel flow using the Launder, Reece, and Rodi 
second-order closure, it must be noted that only mild rotations with Rossby 
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numbers Ro ~ 10 were considered (the reader should note that the Rossby 
number referred to in Launder, Tselepidakis , and Younis 1987 is actually the 
inverse Rossby number). Had those authors considered more rapid rotations 
(i.e., Rossby numbers Ro^O.l) serious inconsistences would most likely 
have arisen as discussed by Speziale (1985). 

Recently, a modified transport equation for the turbulence dissipation 
was proposed by Bardina, Ferziger, and Rogallo (1985) with the purpose of 
accounting for the reduction in dissipation that occurs when isotropic turbu- 
lence is subjected to a rigid body rotation. This model is of the form 
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which differs from the more commonly used model (49) by the addition of the 
last term on the right-hand-side of (56). For isotropic turbulence in a 
rotating frame, Equation (56) takes the form 


de 

dt 


C 4 K C 5 ne * 


(57) 


Bardina, Ferziger, and Rogallo (1985) found that (for = 0.15) Equation 
(57) predicted reductions in the dissipation rate that were in fairly good 
agreement with the experiments of Wigeland and Nagib (1978) for rotating iso- 
tropic turbulence. However, several criticisms can be leveled at Equation 
(56) when applied to anisotropic and inhomogeneous turbulent flows. For 
example, the dependence on the intrinsic mean vorticity term W, W 


does 
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not vanish in a two-dimensional turbulence, thus, violating material frame- 
indifference in this limit. Furthermore, Equation (56) was obtained from (57) 
by an extrapolation which is not unique. More specifically, there are other 
invariants besides ^2 W ij W ij^ which reduce to ft in a rotating 
isotropic turbulence (see Speziale, Mansour, and Rogallo 1987). These prob- 
lems will be addressed in the next section where a hierarchy of consistent 
models will be developed. 
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4. IMPROVED TURBULENCE MODELS FOR NON-INERTIAL REFERENCE FRAMES 

In this section, improved two-equation turbulence models and second-order 
closure models that are consistent with the non-inertial constraints derived 
in Section 2 will be developed. Since rotations can dramatically enhance ani- 
sotropic effects and alter the dissipative properties of the turbulence, eddy 
viscosity models are more likely to yield inaccurate predictions in rotating 
reference frames. Hence, it is best to base two-equation turbulence models on 
some suitable nonlinear generalization of the eddy viscosity models when 
applications to rotating flows are envisioned. Recently, the author developed 
a nonlinear K -fc and K-e model along these lines (see Speziale 1987) 
which appears to account for anisotropic effects much more accurately. This 
model is of the form 


ij 


T K5 iJ 


+ k 1/2 jis 


ij 
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(58) 


in an arbitrary non-inertial frame where C Q is an empirical constant which 
was found to assume an approximate value of 1.68 by correlating with turbulent 
channel flow data. Here, l is the length scale of turbulence which is 
given by 

*3/2 

l - 2C - (59) 

ye 

for the K-e type model. This model constitutes a substantially simplified 
version of a nonlinear eddy viscosity model recently derived by Yoshizawa 
(1984) using Kraichnan's DIA formalism; the simplification primarily arises 
from invoking the constraint of material frame-indifference in the limit of 
two-dimensional turbulence which Yoshizawa's full nonlinear model violates. 
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Since (58) represents a quadratic extension of the linear eddy viscosity 
models which are algebraic in nature, it follows that the satisfaction of 
frame-indifference in the two-dimensional limit restricts any frame-dependence 
in three-dimensional turbulence to arise exclusively from changes in the 
scalar length scale Jt. Unfortunately, such a weak frame-dependence cannot 
account for the considerable additional anisotropies that are caused by a 
moderate system rotation of turbulent shear flows of scientific and engineer- 
ing interest. The constraint of material frame-indifference in the two- 
dimensional limit becomes important in rapidly rotating frames where there is 
a Taylor-Proud man reorganization of the flow — a limit which is largely outside 
of the domain of applicability of such simplified algebraic models which can- 
not account for extremely large anisotropies. Hence, we will relax this con- 
straint in favor of another approximation that follows from a simpified 
analysis of the Reynolds stress transport equation. Bardina, Ferziger, and 
Reynolds (1983) showed, for homogeneous turbulent flows, that the unmodeled 
Reynolds stress transport equation yielded the following analogy: the appli- 
cation of a mean strain in a rotating frame is the same as the appli- 
cation of a mean strain S and mean rotation 28 in an inertial frame of 
reference. This analogy (which is not a rigorous consequence of the Navier- 
Stokes equations since rotational effects arising from the higher-order 
moments were neglected) was shown by Bardina, Ferziger, and Reynolds (1983) to 
be a relatively good approximation for certain rotating turbulent shear flows 
and to be consistent with invariance under the Richardson number. The appli- 
cation of this analogy to the derivations in Speziale (1987) yields a non- 


linear K-e model of the form 
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T iJ - - T K5 iJ + k 1 / Sj - V 2 «ij + s iAj 


(60) 
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is the frame-indifferent Jaumann derivative and the length scale is given by 
i - 2C^K ' /e. For a homogeneous turbulence with constant mean velocity 
gradients in an inertial frame or for turbulent channel flow (the two cases 
considered when the nonlinear K-e model was first calibrated), both (58) 
and (60) reduce to the same form 


T u--! K 5 i j + Kl/ 2 iS i J -V 2 ' s iA j 

(62) 

- -L S S 6,, + to . 1 S, . + w M S, .] 

3 mn mn ij ik kj jk ki 

and hence the value of = 1,68 will not be altered.* It will now be shown 

that this new nonlinear K-e model yields dramatically improved predictions 

for homogeneous turbulent shear flow in a rotating framework (see Figure 4). 

Here, the constant C was taken to be 0.055 (the value recommended by 

y 

Rodi 1984 for homogeneous turbulence where the ratio of the production to dis- 
sipation is equal to two) and the traditional transport Equation (49) for e 


JL 

It is interesting to note that Equation (62) bears a resemblance to the 
nonlinear two-equation models of Pope (1975) and Saffman (1977). 
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was used with C 3 = 1.45 and C 4 = 1.90. A closed form equilibrium solution 
can be obtained which is of the form 


b n ' S c „ 2 # It ' 8 <§»• b 22 ’ V 2 , # t- 4 + 8 <§» < 63 > 
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where the ratio of production to dissipation a = (C^ - 1)/(C^ - 1) and 

2 

the anisotropy tensor = -(t^ + K 6 ) /K • A comparison of the results 

obtained from the linear and nonlinear K-e models (along with the experi- 
ments of Tavoularis and Corrsin (1981) and the large-eddy simulations of 
Bardina, Ferziger, and Reynolds (1983)) are shown in Table 1. Here, the equi- 
librium values of the anisotropy tensor obtained from the nonlinear K-e 
model are dramatically improved with respect to its normal components (the 
reader should note that b^ is not shown since it is precisely equal to 
-(bjj + b 2 2 ) due to the fact that b ^ is traceless). Unfortunately, no 
experimental data is available for rotating shear flow and the values of the 
anisotropies obtained from the large-eddy simulations are somewhat inaccurate 
due to course resolution and the lack of a good defiltering scheme. However, 
there is no question that the normal components of the anisotropy tensor pre- 
dicted by the nonlinear K-e model constitute a considerable improvement 
over their linear counterparts. Both the linear and nonlinear K-e models 
yield the same predictions for the equilibrium values of bj 2 anc * SK/e 
which are the same for all values of ft/S. This is not consistent with 


physical and numerical experiments which indicate that b ^ 2 and SK/e can 
vary considerably with ft/S. In order to predict this dependence, a modi- 
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fied dissipation rate equation must be developed which properly accounts for 
rotational strains — a task of considerable difficulty. 

Now, a consistent modification of the modeled dissipation rate equation 
will be developed which can account for the considerable reduction in dissi- 
pation which occurs in a rapidly rotating isotropic turbulence. As mentioned 

1 i/2 

in Section 3, the intrinsic mean vorticity invariant %^ij^ij^ does 

not vanish in the limit of two-dimensional turbulence. However, this invari- 
ant was arrived at by Bardina, Ferziger, and Rogallo (1985) since it reduces 
to Si in a rotating isotropic turbulence (it was Equation (57) that was used 
to correlate with the experiments of Wigeland and Nagib 1978). Alternatively, 
there are other invariants that reduce to ft for rotating isotropic turbu- 
lence but vanish in the limit of two-dimensional turbulence. The primary such 
invariant is 

I<TB) " It nr + VV 1 1/2 (65) 


which was first introduced by Speziale (1985). This gives rise to the alter- 
native modeled transport equation for the dissipation rate 
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where f^ and f 2 are sufficiently smooth functions of the dimensionless 
invariant I^ Tw ^K/e. For plane homogeneous turbulence, the invariant I^ TW ^ 


reduces to 
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* 2 

i and, hence,for small anisotropies (where K) 

I<TW> : l"i 2 l - <5 Vu )1/2 - (68) 

(69) 

(70) 

(where y^ and y^ are dimensionless constants) it follows that (66) is 
in close approximate agreement with the most recent Bardina modification of 
the dissipation rate transport equation given that y^ = 0.01 and 
y^ * 0.079 (this model also reduces to the more simple model given by (57) 
for rotating isotropic turbulence). Such a model has been shown by Bardina to 
work reasonably well for rotating isotropic turbulence (at moderate rotation 
rates) and for simple plane turbulent shear flows subjected to mild rotational 
strains. However, unlike the Bardina model, this new model for the dissipa- 
tion rate satisfies material frame-indifference in the limit of two- 
dimensional turbulence (as can be seen from (67) since I^ TW ^ + o as 
t 33 * ^ and a H° ws for more general nonlinear dependence on W (the 

simple linear expressions (69)-(70) break down when a wider variety of flows 
is considered). Nonlinear generalizations of (69) and (70) should be pursued 
in future studies. 

Finally, the implications that the non-inertial constraints derived in 
Section 2 have on second-order closure models will be examined. As alluded to 
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before, the Launder, Reece, and Rodi (1975) model as well as the other common- 
ly used second-order closures violate material frame-indifference in the limit 
of two-dimensional turbulence and do not give rise to a Taylor-Proudman re- 
organization for statistically steady turbulent flows in a rapidly rotating 
framework. In fact, Speziale (1985) recently showed that for rotating turbu- 
lent channel flow (see Figure 5) these second-order closure models yield the 
spurious result of a vanishing Reynolds shear stress 
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in the limit as SI +■ « and do not give rise to a full Taylor-Proudman 
reorganization to a two-dimensional state. These problems were not en- 
countered in the recent study of Launder, Tselepidakis , and Younis (1987) 
since they restricted their attention to flows with Rossby numbers greater 
than 10 (a Taylor-Proudman reorganization would only be expected for Rossby 
numbers less than 0.1 — a value nearly two orders of magnitude smaller than 
those considered therein). Complete consistency with the non-inertial con- 
straints (i)-(iv) derived herein can be obtained from second-order closures of 
the general form 
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where e is obtained from the new modeled transport Equation (66) and 

a, 6, y are dimensionless functions of I^ TW ^ and the invariants of r 

(which can be taken to be constants in the first approximation). The first 
term on the right- hand-side of (72), with the coefficient a, arises from 
the rotationally dependent part of the rapid pressure-strain correlation. 
This term was recently derived from a Langevin model by Haworth and Pope 

(1986) who showed that it vanishes in the limit of two-dimensional 

turbulence. The second term on the right-hand-side of (72), with the 

coefficient y» represents the rotationally dependent part of the return 
term which, in a rapidly rotating frame, was shown by Speziale (1985) to give 
rise to a Taylor-Proudman reorganization to a two-dimensional state wherein it 
then vanishes. Here again, * s t1ie third-order diffusion correlation 

which is frame-indifferent and ir^ accounts for the rotationally inde- 

pendent parts of the production, pressure-strain and dissipation rate correla- 
tions. In addition to satisfying material frame-indifference in the limit of 
two-dimensional turbulence (along with consistency with the Taylor-Proudman 
Theorem), this new second-order closure also satisfies constraint (iv). To be 
specific, the rotationally dependent terms in (72) vanish in an isotropic 
turbulence and the modified dissipation rate equation (66) gives rise to re- 
duced dissipation in a rotating frame consistent with constraint (iv). As a 
result of the dramatically improved non-inertial properties of (72), spurious 
physical effects such as (71) (that are predicted by the commonly used second- 
order closures) can be avoided. Equation (72) represents a hierarchy of 
second-order closure models whose detailed study represents an extensive 


research effort that is beyond the scope of the present paper 
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5. CONCLUSION 

In this paper, several important constraints that turbulence models must 
satisfy in non-inertial frames of reference were derived as a rigorous conse- 
quence of the Navier-Stokes equations • Of particular importance was the con- 
straint that turbulence models should only depend on the frame of reference 
through the intrinsic mean vorticity tensor and that all such frame-dependent 
effects must vanish in the limit of two-dimensional turbulence* In addition, 
it was also shown that Rapid Distortion Theory for an isotropic turbulence 
suddenly subjected to a strong rotation can serve as an equally important con- 
straint requiring an initially isotropic turbulence to decay isotropically 
(with a reduced dissipation rate) in a rotating frame* All of the commonly 
used turbulence models were shown to be in serious violation of these con- 
straints and, thus, inconsistent with the Navier-Stokes equations. An im- 
proved two-equation turbulence model was developed which was demonstrated to 
be substantially superior to the more standardly used K-e model in the 
description of homogeneous turbulent shear flow in a rotating frame. Further- 
more, a hierarchy of consistent second-order closure models was developed 
which have dramatically improved properties in rotating frames over the more 
commonly used second-order closures* A complete calibration and testing of 
such models is a massive research effort that is beyond the scope of the 
present study. However, such work is currently underway in collaboration with 
others. 

Finally, it should be mentioned that the results of this study could have 
important implications in the analysis of curved turbulent flows. As demon- 
strated herein, once the inertial form of a turbulence model is specified, its 
non-inertial form is automatically determined by appropriately replacing the 
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mean vorticity with the intrinsic mean vorticity. Consequently, if a turbu- 
lence model exhibits incorrect behavior in a non-inertial frame, this means 
that the dependence of the inertial form of the model on the mean vorticity is 
faulty. Since the mean vorticity plays an important role in the description 
of curved turbulent flows, it is quite likely that the difficulty in describing 
such flows is a result of the use of models that exhibit physically incorrect 
non-inertial behavior. A more detailed discussion of this point will be the 
subject of a future paper. 
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■ 

Equilibrium 

Linear 

Nonlinear 


Large-Eddy 

1 

Values 

K-e Model 

K-e Model 

Experiments 

Simulations 


b ll 

0 
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0.403 

0.606 ! 

I 

i 

o 

II 

b 22 

0 

-0.308 

-0.295 

1 

-0.530 ! 
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1 

CO 

c: 

b 12 

-0.332 

-0.332 

-0.284 

-0.288 

i 


SK 

e 

6.03 

6.03 

6.08 



b ll 
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0.0616 

- 

0.121 

in 

CM 

b 22 

0 

0.0616 

- 

0.091 

1 
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J 

It 

CO 

C! 

b 12 

-0.332 

-0.332 

- 

-0.697 


SK 

e 

6.03 

6.03 

- 

- 


b n 

0 

-0.308 

- 

-0.530 

0.5 

b 22 

0 

0.431 

- 

0.500 

H 






CO 

c: 

b 12 

-0.332 

-0.332 

- 

-0.197 


SK 

e 

6.03 

6.03 

- 

- 


Table 1. Equilibrium results for homogeneous turbulent shear flow in a ro- 
tating frame: Comparison of the predictions of the K-e model with the 
large eddy simulations of Bardina, Ferziger, and Reynolds (1983) and the ex- 
periments of Tavoularis and Corrsin (1981). 
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Figure 5. Turbulent channel flow in a rotating frame 
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